Current and spin torque in double tunnel barrier 
ferromagnet— superconductor— ferromagnet systems 
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We calculate the current and the spin torque in small symmetric double tunnel barrier 
ferromagnet-superconductor-ferromagnet (F-S-F) systems. Spin accumulation on the supercon- 
ductor governs the transport properties when the spin-flip relaxation time is longer than the trans- 
port dwell time. In the elastic transport regime, it is demonstrated that the relative change in the 
current (spin torque) for F-S-F systems equals the relative change in the current (spin torque) for 
ferromagnet-normal metal-ferromagnet (F-N-F) systems upon changing the relative magnetization 
direction of the two ferromagnets. This differs from the results in the inelastic transport regime 
where spin accumulation suppresses the superconducting gap and dramatically changes the magne- 
toresistance [S. Takahashi, If. Imamura, and S. Maekawa, Phys. Rev. Lett. 82, 3911 (1999)]. The 
experimental relevance of the elastic and inelastic transport regimes, respectively, as well as the 
reasons for the change in the transport properties are discussed. 

PACS numbers: 72.25.Mk, 74.50.-|-r, 75.70.Pa, 73.40.Gk 



I. INTRODUCTION 

Conventional low-temperature superconductors result 
from the pairing of electrons with spin up and spin down 
and opposite momentum. Ferromagnets in contact with 
a nonmagnetic metal can induce a nonequilibrium spin 
accumulation in the nonmagnetic metal. Such a spin 
accumulation in a superconductor can reduce the acces- 
sible number of pairs of spin-up and spin-down electrons. 
Consequently, spin accumulation in superconductors can 
reduce the superconducting gap and dramatically change 
the transport properties. 

Spin accumuLation and its effect oii,the current-voltage 
characteristicalfa and the shot noiseEl have been studied 
thoroughly in double barrier F-N-F systems. We will 
here consider double tunnel barrier F-S-F systems with 
general noncollinear magnetization directions. The influ- 
ence of spin accumulation on the superconducting gap, 
the current and the spin torque will be considered. 

Spin accumulation and its influence on the supercon- 
ducting gap strongly depend on the competition between 
different relevant transport time scales; the transport 
dwell time characterizes the typical time an electron 
spends on the superconducting island on passing through 
the system, the spin-flip relaxation time Tgt gives the 
time-scale for the coherence of the electron spin, and the 
energy relaxation time te denotes the time scale for the 
interchange of energy between the quasiparticle and the 
rest of the system. Spin accumulation in F-N-F or F-S- 
F systems requires that the spin-flip relaxation time is 
much longer than the transport dwell time, Tgt r^. The 
transport dwell time decreases with decreasing size of 
the system via the density of states and consequently 
spin accumulation will only take place in mesoscopically 
small systems (for a discussion, see Ref. |^). F-SiF sys- 
tems have been studied in an interesting articleQ in the 
limit when Tgf ^ Td ^ te, i.e., when the energy relax- 
ation is sufficiently strong so that the quasiparticles for 



each spin relax to Fermi-Dirac (FD) distribution func- 
tions with spin-dependent chemical potentials. Wc will 
below show that a different regime {TE,Tsf) 3> Td is rel- 
evant in some experimental situations, and we will work 
out the transport properties and compare them to the 
inelastic regime. 

Our aim is to provide a complete description of the 
elastic transport regime for F-S-F systems with non- 
collinear magnetization directions. We will alsa,general- 
ize the results in the inelastic transport regimeu to non- 
collinear magnetization configurations. Both the current 
and the spin torque relevant for magnetic switching be- 
havior in small ferromagnetic particlesErB will be studied. 



II. TRANSPORT REGIMES 

The relevance of the elastic or inelastic transport 
regimes depends on two time scales intrinsic to the su- 
perconducting island, the energy relaxation time te and 
the spin-flip relaxation time Tgt as well as the time-scale 
governed by the contacts to the superconducting island, 
Td. As discussed in the Introduction, spin accumulation 
requires Tgt ^ Td and we will assume that this is satis- 
fied in the following. This implies that when the intrinsic 
time-scales are in the regime te ^ Tgf, the spin accumu- 
lation requirement dictates the elastic transport regime, 
Te 3> Td. In the opposite limit, te Tgf both elastic 
and inelastic regimes can be relevant depending on the 
contacts. Let us now discuss the ratio between the spin- 
fiip relaxation time Tgf and the energy relaxation time 
Te for conventional s-wave superconductors. First, let us 
consider the case of Al. 

The energy relaxation rate 1 /te has contributions due 
to the electron-electron scattering (1/Te_e) and due to 
the electron-phonon processes, i.e., quasiparticle- phonon 
scattering and quasiparticle recombination with phonon 
emission to form Cooper pairs, (l/re_ph). It is usually 
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assumed that the two contributions add up: 1/te = 
1/Te_e + 1/Te_ph- In most Conventional superconduc- 
tors, the dominant quasiparticle energy relaxation rate is 
due to the electron-phonon processes, except for metals 
with a large Debye frequency and a small superconduct- 
ing gap, such as Al and Zn, for which-|electron-clcctron 
scattering processes can be significant.Ll 

The energy relaxation time due to the electron-phonon 
processes was studied in Ref. ^ for a variety of materials. 
For Al, Te_ph « 5 X 10~* s close to the Fermi level at the 
critical temperature. A further reduction of the tempera- 
ture leads to an even smaller rate for the electron-phonon 
processes at the relevant energy A(r)E Furthermore, the 
electron recombination rate to form Cooper pairs is re- 
duced for spin-polarized quasiparticles. 

The energy relaxation time due to the electron-electron 
interaction in dirty normal state [Systems at an energy e 
with respect to the Fermi level isO 
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where kp is the Fermi wave vector, r is the elastic scatter- 
ing time, and A = vpT is the mean free path. Relevant 
energy scales for F-S-F systems are around the super- 
conducting gap (in Al, « 1.2 K) and lower. Using a 
typical Fermi energy for metallic systems Ep = 10 cV 
and a value of the electron mean free path for a dirty 
system kpX = 10 (which corresponds to A ~ 6 A) we 
find for Al r^-e ~ 4 x 10^^ s. For cleaner systems, re_e 
can be much larger since Te_e cx: r^^^. 

The spin-flip relaxation time depends on the relativis- 
tic spin-orbit interaction and the magnetic impurity scat- 
tering. For nonmagnetic impurities, the spin-orbit in- 
teraction dominates the spin-flip rate. The spin-orbit 
scattering time in superconducting Al was found to be 
Tso « 2x 10~^^ s in the pioneering tunneling spectroscopy 
measurements by Tedrow and Meservey.a In their experi- 
ment they used an Al film with the elastic scattering time 
T K, 10^^*^ s (corresponding to A w 20 nm) from which we 
can estimate Te_e using Eq. (|l|). For such a clean system, 
the electron-electrorLinteraction is weak and te ~ Te-ph. 
Johnson and SilsbeeEEl subsequently found a longer spin- 
flip relaxation time in an Al (normal state) thin film, 
Tst ~ 10~^ s. To the best of our knowledge there are no 
other experiments that have been so successful in achiev- 
ing such a long spin-fiip relaxation time in Al. Besides, 
the spin-flip relaxation time decreases in small metal paii-i 
tides due to the enhanced scattering rate at the surfacecl 
and in materials wiiii magnetic impurities due to the 
spin-spin scattering.tll We, therefore, consider the long 
spin-flip relaxation Tgf ~ 10^* s an optimistic estimate 
in small superconducting Al particles and believe that 
times of the order rgf ~ 10~^^ s are currently easier to 
achieve experimentally. 

The energy relaxation time can therefore be of the 
same order as the spin-flip relaxation time or larger for 
superconducting Al particles. It is therefore likely that at 
least some experiments on F-S-F systems with Al island 



will be in the regime te ^ Tgf ^ r^. This motivates our 
study of the elastic transport regime. 

Let us extend this discussion to the case of some other 
common s-wave superconductors. Kaplan et alB showed 
that for many simple conventional superconductors (Al, 
Zn, Nb, In, Sn, Ta, Hg, Tl, Pb) Te^ph ranges between 
10"^ s and 10~^^ s. We have to compare it to the spin- 
orbit scattering time Tso- Approximately, Tso scalesJin- 
early with the elastic scattering time r according tdlj 



(2) 



where a is the fine structure constant and Z is the atonuc. 
number of the impurity atoms. Meservey and Tedrowlij 
have showed by compilation of data that Eq. (||) is a rea- 
sonable estimate of the spin-orbit scattering rate for very 
clean small particles or films, when the surface collisions 
dominate the elastic scattering processes. Nevertheless, 
they showed that for many metals with atomic numbers 
under 100, Eq. (||) typically gives an underestimate of 
the scattering rate 1/tso by a factor of the order of 10. 
According to Eq. (||), in metals heavier than Al, it is even 
more difficult to achieve a lower spin-flip relaxation rate. 

We therefore conclude that the elastic regime 
{TE,Tsf) 3> Td can be achieved in experiments on F-S-F 
systems with a variety of materials used for the supercon- 
ducting node. This is true for relatively clean systems, 
when Te-ph «C Te-e- FoT dirty systems, electron-electron 
interactions can become significant and since Tf.-e oc r^/^ 
(|l|) and Tso r (||), Tp-e/Tso ^/t in the theoretical 
limit r — > 0. If this limit is physically achieved, both elas- 
tic and inelastic regimes can be relevant. The above dis- 
cussion does not apply to high- Tc superconductors which 
typically have much higher energy xelaxation rate. In 
particular, some recent experiments^ on spin-polarized 
quasiparticle injection in the cuprate superconductors in- 
dicate the relevance of the inelastic regime as studied in 
Refs. I and 
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In the elastic regime, the ratio between the energy- 
relaxation time and the spin-flip relaxation time does 
not matter; the occupation of the quasiparticles states 
are not determined by FD distributions, but has to be 
determined from the transport equations. These trans- 
port equations will be found in the following. 



III. MODEL 

The F-S-F system consists of a superconducting island 
connected to a left and a right ferromagnet by symmetric 
tunnel junctions. The left and the right ferromagnets are 
attached to a left and a right reservoir by good metallic 
contacts; it is assumed that the current is only limited 
by the tunnel conductance between the ferromagnets and 
the superconductor. Each monodomain ferromagnet has 
a well-defined magnetization direction, m/ and ra.r for the 
left and the right ferromagnet, respectively. The system 
is biased by an external voltage source between the left 
and the right reservoirs. The convention is such that the 
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chemical potential in the left reservoir is higher than that 
in the right reservoir by eV (e > 0). We will consider 
the current and spin current from the superconducting 
island to the left reservoir, and the spin torque on the 
left ferromagnet. 

We consider the regime when the tunnel conductance 
is much larger than the quantum conductance so that 
the Coulomb charging effects can be disregarded. It is 
further assumed that the level spacing is much smaller 
than the temperature and the superconducting gap, so 
that the superconducting state can be well described by 
the BCS theory. The dwell time is much larger than the 
time for the diffusion through the superconducting island 
so that the phase-space occupation of the quasiparticles 
can be described semiclassically. It is also assumed that 
the spin-flip relaxation time is longer than the transport 
dwell time in all calculations, Tgf S> t^, so that a nonequi- 
librium spin accumulation on the superconducting island 
can exist. 

We describe the electron transport through a tunnel 
junction sandwiched between a superconductor and a 
ferromagnet using a phenomenological tunneling Hamil- 
tonian. The total Hamiltonian of the system is H — 
Hq + H', where 



where 



^0 = X! ^i^kCLk + 



(3) 



is the Hamiltonian of the uncoupled superconductor and 
ferromagnet and 



H' = Y^ 4^kpak + H.c. 

kp 



(4) 



is the phenomenological tunneling Hamiltonian. The fer- 
romagnet has a general magnetization direction m. Cp 
(ofe) is a column vector of spin-up and spin-down annihi- 
lation operators for an electron with a momentum p (k) in 
the ferromagnet (superconductor), T^p = Tkpiii^ +Tkpiu^ 
and ip — EpiM^ -|- CpiiiK The tunneling matrix T and the 
ferromagnet energy matrix e depend on the magnetiza- 
tion direction of the ferromagnet through the projection 
matrices ii^ = (l + & ■ m) /2 and = (l — & ■ m) /2, 
where m is a unit vector in the direction of the local 
magnetization and & = {d'x,(^y,^z) is a vector of the 
2x2 Pauli matrices. 

The 2x2 current matrix in the spin space evaluated 
inside the superconductor is 



fa/3 



= e-Tr(AfT)+H.c. , (5) 



where M^^^ = 



and the trace is taken over the 
momentum space indices. The current is / = 



and the spin current is Is = —h/ (2e) X^a^ ^ 



'af3 



In order to calculate Af , we express the electron oper- 
ators a in terms of the quasiparticle operators 7 by the 
Bogoliubov transformation: 



1 
1 



Vk = Vk 



1 

-1 



and Uk = [{l+^k/Ek)/2f\ Vk = [{1 - ^k/Ek)/2f\ 
where is the electron energy relative to the chemical 
potential, Ek — + is the corresponding quasi- 
particle energy and A is the superconducting gap. 

Up to this point, everything is exact within the tunnel- 
ing Hamiltonian phenomenology. We calculate M [and 
thus find the current (|^)] to the lowest nonvanishing or- 
der in T. This approximation neglects Andreev reflection 
and higher-order corrections to the quasiparticle current. 
In the next section, we show the final result for the cur- 
rent (|^) and apply it to describe the transport properties 
of double tunnel barrier F-S-F systems. 



IV. THEORETICAL RESULTS 

The current operator is a 2 x 2 matrix in the spin 
space. Let us first include the result when the islaiid. 
is in the normal state. It was recently demonstratedll3 
that the 2x2 current (per unit of energy, at a given 
energy e) through a single normal metal-ferromagnet (A^- 
F) junction at the normal metal side can be written as 

IN (e) e = g^ii'^ S/n (e) + g^u^ SJn (e) 

+g^^v) Sfr, (e) + g^^u^ Sf^ (e) v) , (7) 

where 5fN (e) is the 2x2 matrix difference in the dis- 
tribution function on the ferromagnetic and the normal 
metal sides 5fM (e) = fp (e) — In (e)- The semiclassical 
distribution function / (e) is defined by 



(aL.(e)a/3n {e')) = &nrn5 {e - e') [e) 



(8) 



Ofe = Uklk + Vk"f^k ; 



(6) 



where aj^^ creates a spin-a particle in the m"" quan- 
tum state. In Eq. (|^), g^ {g^) denotes the junction 
conductance for a spin-up (spin-down) electron, and g^^ 
(= [g^^Y) is the mixing conductance first introduced in 
Ref. |l6[ For a tunnel junction, as we study here, we find 
that the calculation based on the tunneling Hamiltonian 
(^) gives the mixing conductance g^^ = (g^ -|- g^) /2, 
which is consistent with the alternative tunneling Hamil- 
tonian phenomenology in Ref. 

The current through the F-N-F system is found by us- 
ing the charge and spin conservation on the normal metal 
node, so that the current through the first junction can- 
cels the current through the second junction; in the elas- 
tic transport regime, the charge and spin arc conserved 
at each energy: «i(e)-l-i2(e) = 0, and in the inelastic scat- 
tering regime, the distribution functions are assumed to 
be FD distributions for each spin direction and the total 
charge and spin are conserved: de [ii(e) -I- 12(e)] = 0. 
By using the current conservation, the distribution on the 
normal metal node can be found, and, by inserting the 
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resulting distribution into Eq.- 
through the tunnel junction isE 



n), the resulting current 



I{V) 



1 



p2 sin^ 61/2] 



5T +5i 



V . 



(9) 



The current only depends on the total conductance of 



spin- up and spin-down electrons + and the relative 
polarization of the contacts p = {g^ — 5x)/(ffT 9i) ^'^'^ 
has a simple |Sin^ 9/2 depends on the relative magnetiza- 
tion angle 6'.E3 For double barrier F-N-F systems with 
general contacts a macfi general and complicated angular 
dependence is found.E^I 



J 



In this paper we generalize Eq. (0) to describe quasiparticle current through conventional s-wave superconductor- 
ferromagnet {S-F) junctions. We find that the quasiparticle current iqp can similarly be expressed by Eq. (|^) by 
replacing S/n with an effective difference in the distribution functions 6fs and renormalizing the conductance by the 
BCS density of states: 



V(e)e = Ns{e) g^u^6fs{e)u'^ + g^u^Sfs{e)u^ + g'^^u'^ Sfs{e)u^ + g^'^u^Sfs{e)u'^ 



where 



Sfs (e) 



(6) - 1 + i ( T> i-e) + T< (-6) j + ^ (f> i~e) - T< {-e) 



e> A 
e < -A 



(10) 



(11) 



Here, A^s(e) = jej/Ve^ — A^ is the normalized BCS density of states, Q(e) = Ve^ — A^/e is the quasiparticle effec- 
tive charge, fg (fg) is the superconducting distribution function of electron-like (hole- like) quasiparticles (in the 
superconductor, the distribution function is defined via quasiparticle operators jam related to electron operators aam 
by the Bogoliubov transformation (^)). is a "flipped" distribution function /; that is if / = fol + /s<t • u, then 



/qI — /sir • u. The energy e is measured with respect t o the chemical potential of the Cooper pair condensate. 



It is important to note that e in Eqs. (|T^) and ( [III ) is 
the electron energy and thus can be positive or negative; 
the quasiparticle energy in the superconductor (which we 
will denote E), on the other hand, is always positive. 

The superconducting gap is determined by the 
(nonequilibrium) occupation of the quasiparticles: 



Tr fs{E) 
dE^^= 

A VE^ - A- 



(12) 



where fs{E) = fs{E) + fsi^) /2 is the distribution 
function for combined electron-like and hole-like excita- 
tions and Aq is the zero-temperature equilibrium gap. 

In the elastic case {te ^ Td), we solve for the distribu- 
tion function in the superconductor based on the conser- 
vation of charge and spin on the superconducting node 
at every energy e. The resulting distribution function 
depends on the spin-up and spin-down junction conduc- 
tances, voltage bias, superconducting gap A^i and the 
relative magnetization direction of the ferromagnets. On 
the other hand, we find that the average occupation num- 
ber of a particular quasiparticle energy level E > Agi on 
the superconducting island is simply 

T4fsAE)] = f{E - eV/2) + f{E + eV/2) , (13) 

where f{E) is the Fermi function. Eq. ( p^ ) is a direct 
consequence of the charge conservation at electron ener- 
gies e = ±E. 



The gap Ad is found self-consistently by inserting 
Eq. (P) into Eq. (Ill 



In 



Ao 
Aci 



dE 



X [f{E~eV/2) + f{E + eV/2)] .(14) 

The gap is thus independent on the magnetization di- 
rections and the polarization of the junctions. This is 
very different from the results in the inelastic trans- 
port regime studied in Ref. ^. The gap equation ( |l^ ) 
gives the same gap suppression as for normal metal- 
superconductor-normal metal (N-S-N) systems, e.g. as 
for a F-S-F system with no spin polarizativeoi. The prop- 
erties of N-S-N systems are well known:0 At T = 0, 
for \eV/2\ < Aq, electrons cannot tunnel onto the su- 
perconducting island and the gap is invariant, but for 
\eV/2\ > Aq, electrons tunnel and induce an excess num- 
ber of quasiparticles on the island which destroy the su- 
perconductivity altogether. Consequently, at T = 0, 
Aci(V^) = Aoe(Ao - \eV/2\), where 9 is the Heavi- 
side step-function. Increasing the temperature decreases 
Aoi(V^) and the critical voltage for the superconductor- 
normal metal phase transition. Moreover, Aci(T^) re- 
mains single- valued. This is in contrast to the inelastic 
transport regime when the gap Ain(T^) is multiple- valued 
for some values of V and T and the system exhibits hys- 
teresis, so that the thermodynamic properties of the is- 
land depend not only on the applied voltage but also on 
the history of the system. 
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In the elastic transport regime, the magnetization di- 
rection dependent quasiparticle current is 

/ei(0, V, T) = [1 - p2 sin2 e/2] /o {V, T, Aei(F, T)) , 

(15) 

where /o(V, T, A) is the current in the parallel configura- 
tion [9 = 0) of the magnetization directions and the gap 
Aoi has to be determined self-consistently at each voltage 
bias according to Eq. ([ij) . In the parallel magnetization 
alignment, there is no spin accumulation and the current 
/o is given by the familiar result for N-S-N systems with 
the total N-S junction conductance + g^: 

/o(l/,T,A) = {g^+g^)/e / dENs{E) 

J A 

X [f{E - eVl2) - f{E + eV/2)] . 

(16) 

Let us compare the results above in the elastic trans- 
port regime to the transport properties in the inelas- 
tic regime [te <C t^). The inelastic transport regime 
has been|-studied when the magnetization directions are 
collinearB It was shown that in the parallel (P) align- 
ment of the magnetizations, there is no bias dependence 
of the gap (that is the gap Ain always equals its equilib- 
rium value A(T)) and the quasiparticles are distributed 
according to the equilibrium FD distribution. In the an- 
tiparallel (A) alignment, on the other hand, the gap is 
suppressed.cl This is because the total (integrated over 
all energies) spin and charge conservation on the super- 
conducting island requires the chemical potentials of the 
spin-up and spin-down quasiparticles to be shifted oppo- 
sitely by a finite value with respect to the chemical 
potential of the condensate. S^i/Aq is a function of V/ Aq, 
T/Aq and the ferromagnet polarization p. Knowledge of 
this function gives the full thermodynamic and transport 
description of the F-S-F system in this configuration.a In 
particular, 

Tr[/s,in(i?)] = fiE - 5ii) + f{E + 6ii) , (17) 

in contrast to Eq. (|l3|), and the gap Ain(V, T, p) can be 
found from the self-consistency equation (|l2|). 

Let us generalize the result in Ref. ^ to noncoUinear 
magnetization configurations. For an intermediate align- 
ment of the ferromagnet magnetizations, we can find the 
magnitude of spin accumulation dfi and a unit vector 
u (the direction of the nonequilibrium magnetization) 
which diagonalizes the quasiparticle distribution function 
on the island in the spin space, by requiring conservation 
of the total spin and charge on the superconducting node. 
We find that u points along the relative magnetization 
of the two ferromagnets, u cx (m; — m^) and Sfi (and, 
therefore, Ajn) depends on the polarization p and rela- 
tive angle 6 through the combination psin6'/2. In other 
words, increasing angle with respect to the A alignment 
at a fixed bias and temperature has the same effect on 
the thermodynamic properties of the island as decreasing 
the polarization while keeping the alignment fixed. 



The result for the quasiparticle current can than be 
written similarly to the elastic case ( |l5| ) as 

Iin{0,V,T) = [1 - p2 sin2 (?/2] 

x/o(V^,T,Ai„(F,r,psin^/2)) . (18) 

The functional form of the gap Ain{V, T,p) has been de- 
scribed in detail in Ref. || for the antiparallel alignment 
and this can be used for the general situation on let- 
ting p — > psin6'/2. Comparing Eqs. ( p^ and (|lj), one 
can see that the only difference between the elastic and 
inelastic regimes in this context is the different supercon- 
ducting gap Aoi and Ajn, respectively. While in the elas- 
tic regime, gap suppression is independent of the relative 
angle 0, in the inelastic regime, it depends on the angle 
through the combination p sin 9/2. Keeping this distinc- 
tion in mind, we can drop the arguments in Eqs. ( p^ and 
( [Tsl ) and rewrite them as 

/= [l-/sin2 6l/2]/o. (19) 

An important observation to make is that for a given 
junction conductance, magnetization configuration, volt- 
age bias and temperature, the gap suppression is greater 
in the elastic than in the inelastic regimes. This follows 
from Eq. ( |l3| ) for the elastic case, Eq. (|l^ for the inelas- 
tic case, the self-consistency requirement (|l^ and the 
fact that Sn < eV/2E 

We can summarize the results in the elastic (|l^) and 
inelastic regimes ( p^ by noting that the gap suppression 
and the tunnel conductance of a symmetric F-S-F system 
in the noncoUinear alignment can be found by calculating 
the corresponding values in the A alignment and using 

Pcff=P sin 61/2 (20) 

instead of p. It is true for both elastic and inelastic 
regimes. In particular, Eq. (|20| ) implies that in the in- 
elastic regime, the superconducting island undergoes a 
superconductor-normal metal phase transition on vary- 
ing the relative magnetization direction gradually from 
the P to A alignment at a fixed voltage. This occurs when 
the ferromagnet polarization is large enough so that the 
superconducting gap vanishes in the A configuration. 

Spin-polarized current driven through a magnetic mul- 
tilayer system can induce torques on the magnetic lay- 
ers due to the spin transfer. Many new phenomena re- 
lated to these current-driven spin torques in magnetic 
layered sysjtems were both predicted and observed in re- 
cent years. Bfl'tS Related to our subject, for example, is 
torque on the ferromagnet layers in five-layer N-F-N-F-N 
systems studied in Ref. ^, where it was shown how the 
spin torque can be responsible for the mesoscopic preces- 
sion and switching of the magnetization directions of the 
ferromagnet layers. 

We calculate the torque on the ferromagnets when the 
spin current passes through the double tunnel barrier F- 
S-F system. The torque exerted on the left ferromagnet 
is 

X = t - (t • mOm, . (21) 
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First, we find the spin current for the F-S-F system 
(by spin conservation on the node, it is the same in both 
junctions): 



/ , \^'P T 
-(nir + m;)— /o- 

4e 



(22) 



Here, Iq is given by Eq. (jLuj), and the gap A has to 
be evaluated self-consistently, as described above for the 
elastic and inelastic transport regimes. The sign conven- 
tion for the current and the spin current in E q. (22) is 
described in the first paragraph of Section 
From Eqs. (|2l|) and ^ it follows that 
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^ hp sin 9 



(23) 



where the unit vector r = (m; cos 9 — nir)/ sin 9. 

In the elastic voltage-biased system, the gap A and, 
consequently, Iq (|l^) do not depend on 0. Therefore, the 
torque t{9) ( p3| ) has a simple sinusoidal dependence on 
the relative angle 9; in particular, t{9) reaches its max- 
imum when 9 — tt/2. However, in the inelastic trans- 
port-regime, the gap suppression is larger with larger 
Pcff The gap is most suppressed in the A configuration 
{PcS = p) and is not suppressed at all (with respect to its 
equilibrium value A(T)) in the P configuration (jjcs = 0). 
On the other hand, according to Eq. (p^), higher gap 
suppression means higher current Iq and, consequently, 
higher torque r oc /q (p3[). As a result, the maximum 
in torque t{9) is shifted to 6* > tt/2 in this regime. In 
general, for a given junction conductance, magnetization 
configuration, voltage bias and temperature, the torque 
will be larger in the elastic than in the inelastic regime 
as the gap suppression is greater in the elastic regime. 

The torque-to-current ratio can be found by using 
Eqs. (§|) and (|l|): 



r/I 



psmt 



4e 1 -p2 sin^ g/2 



(24) 



In the case of semimetallic contacts, 0-,= 1, Eq. (24) is 
the same as the result of SlonczewskiB for N-F-N-F-N 
systems: 



t/I = — tan6'/2, 

even though the angular dependence of the current / is 
different in our model. 



V. CONCLUSIONS AND DISCUSSIONS 

Spin accumulation is crucial for the_superconducting 
gap suppression in the inelastic regime.cl The mechanism 
of this suppression was discussed in detail in Ref. ^ for 
the case of antiparallel alignment and we generalized this 
to noncoUinear magnetization configuration. 

In contrast to this, in the elastic regime, spin accu- 
mulation does not affect the superconducting gap. This 



follows from Eq. (p^: even though the distribution of 
the quasiparticles is anisotropic in the spin space for fi- 
nite ferromagnet polarization p, the polarization p does 
not affect the total distribution of the quasiparticles in 
the momentum space. This is a direct consequence of the 
charge conservation on the superconducting island at ev- 
ery energy level. In this transport regime, the current 
is proportional to the current in a similar N-S-N sys- 
tem and the relative current change upon changing the 
magnetization directions is identical to the situation of a 
normal metal island. 

We derived a very simple result for the spin torque in 
F-S-F systems ( 23 ) which can be viewed as a generaliza- 
tion of the result for F-N-F systems. In the elastic case, 
the relative torque change upon changing the angle 9 is 
identical to the situation of F-N-F systems. In the inelas- 
tic case, the angular dependence of the torque becomes 
more complicated as the gap suppression is angular- 
dependent. The maximum torque sets in for 9 > 7r/2 
due to the suppression of the gap that increases the spin 
current. 

Let us briefly discuss generalizaiians to go beyond the 
tunnel regime. It has been shownta that Eq. (0) is the 
general expression for the current through an arbitrary 
N-F junction when the spin-flip relaxation processes can 
be disregarded. We believe that Eq. (^0|) also holds gener- 
ally for the quasiparticle current iqp in S-F systems with 
arbitrary junctions if we redefine Ng{e) to be an energy- 
dependent normalization factor, which equals normalized 
ECS density of states for the tunnel contacts but has to 
be found separately in the case of other contacts. For 
higher transparency junctions, Andreev reflection also 
participates in the transport properties. The Andreev 
current is spinless since a Cooper pair consists of an elec- 
tron and a hole with opposite spins, 

iA{e)e = igA{e/A) (tv [^(e) + ^(-e)] - 2) , 

where gA is an energy-dependent conductance. The total 
current is given hy i — iqp -\- ia ■ 

In conclusion, we have discussed the relevance of the 
elastic and inelastic transport regimes and showed that at 
least some experiments on F-S-F systems with relatively 
clean s-wave superconducting island can be in the elastic 
regime. 

We generalized the quasiparticle current for N-F tun- 
nel junctions (^ to describe S-F tunnel junctions by 
defining an effective difference in the distribution func- 
tion across the junction (|lll). We then used it to find the 
thermodynamic gap suppression, electric current-voltage 
characteristics and mechanical spin-torque properties of 
symmetric double barrier F-S-F systems biased by an ap- 
plied voltage source. We considered both the elastic and 



re inelastic regime has 
in the case of collinear 



the inelastic transport regimes, 
been described in detail elsewheri 
magnetization directions of the ferromagnetic reservoirs 
and here we generalized it to the case of arbitrary align- 
ment of the ferromagnets. Angular dependence of the su- 
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perconducting gap suppression and current-voltage char- 
acteristics can be succinctly described using effective fer- 
romagnet polarization (pO|). In the inelastic transport 
regime, Eq. (|2^) can be used to determine the angle of 
the phase transition from the normal to superconducting 
state when the alignment is gradually varied from the 
antiparallel to parallel at a fixed voltage bias. This oc- 
curs when the polarization is large enough so that the 
island superconductivity is completely suppressed in the 
antiparallel configuration. 
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